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We study the one-dimensional S = 1 Blume-Emery-Grifliths model. Upon transforming the 
spin model into an equivalent fermionic model, we provide the exact solution within the Green's 
function and equations of motion formalism. We show that the relevant response functions as well as 
thermodynamic quantities can be determined, in the whole parameters space, in terms of a finite set 
of local correlators. Furthermore, considering the case of an antiferromagnetic chain with single-ion 
anisotropy in the presence of an external magnetic field, we find three plateaus in the magnetization 
curve; in the neighborhood of the endpoints of the intermediate plateau, the specific heat shows a 
double peak structure. 
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I. INTRODUCTION 

The spin-1 Ising model with bilinear (J) and bi- 
quadratic {K) nearest-neighbor pair interactions and a 
single-ion potential (A) is known as the Blume-Emery- 
Griffiths (BEG) model [l|. With vanishing biquadratic 
interactions, the model is known as the Blume-Capel 
model I,!- 

There is a large and diffused interest in the study of this 
model, motivated by several reasons. The BEG model 
was originally introduced to describe the phase separa- 
tion and superfluidity in the '^He-'^He mixtures, but it can 
also describe the properties of a variety of systems rang- 
ing from spin-1 magnets to liquid crystal mixtures, mi- 
croemulsions, semiconductor alloys, to quote a few. Both 
the BEG and BC models have been investigated using 
many different approaches for different lattice type and 
dimensions. In one dimension and zero magnetic field, 
the spin-1 Ising model and the BEG model have been 
solved exactly by means of the transfer matrix method 
P] and by means of the Bethe method @- Exact so- 
lutions have also been obtained for a Bethe lattice [y| 
and for the two-dimensional honeycomb lattice Q- For 
higher dimensions, among the various approximate and 
simulation techniques, the most common approach to the 
BC and BEG models is based on the use of mean field 
approximation. However, renormalization group studies 
show some qualitative differences from the mean field re- 
sults. The one-dimensional case for the BEG model was 
studied in Ref . [8| , where exact renormalization-group re- 
cursion relations were derived, exhibiting tricritical and 
critical fixed points. We refer the interested reader to 
Ref. [9] for a rather broad list of works devoted to the 
study of the BC and BEG models. 

In previous works 'l(|, we have shown that, upon 
transforming to fermionic variables, spin systems can be 
conveniently studied by means of quantum field meth- 
ods, namely: Green's functions and equations of motion 



methods. This approach has the advantage of offering a 
general formulation for any dimension and to provide a 
rigorous determination of a complete set of eigenopera- 
tors of the Hamiltonian and, correspondingly, of the set 
of elementary excitations. In this paper we apply this for- 
mulation to the ID BEG model, and we obtain the exact 
solution of the model in the whole space of the parame- 
ters. Postponing to a forthcoming work a comprehensive 
analysis of the model, in this article we focus our study on 
the thermal and magnetic properties of the antiferromag- 
netic (J < 0) spin-1 chain with single-ion anisotropy in 
the presence of an external magnetic field, in the limit of 
vanishing biquadratic interaction. This particular case is 
indeed interesting since one finds magnetization plateaus, 
experimentally observed [l^l- At T = the magnetiza- 
tion curve forms plateaus with abrupt jumps from one to 
another at certain values of the magnetic field h. When 
A = 0, the ground state is purely antiferromagnetic for 
— 2|J| < h < 2\J\. By varying the magnetic field, the 
system undergoes a phase transition to a pure ferromag- 
netic regime at h — ±2|J|. When the anisotropy A is 
turned on, one observes an intermediate phase charac- 
terized by half of the spins oriented along the external 
field and the rest perpendicular to it. The intermediate 
phase between the anti- and ferromagnetic ones, has a 
width depending on A whose endpoints are denoted by 
he and hs, i.e., critical and saturated field, respectively. 

The paper is organized as follows. In Section [TTl upon 
introducing a complete set of composite operators, eigen- 
operators of the Hamiltonian, we outline the analysis 
leading to the algebra closure and to analytical expres- 
sions of the retarded Green's functions (GFs) and corre- 
lation functions (CFs). Since the composite operators do 
not satisfy a canonical algebra, the GF and CF depend 
on a set of internal parameters, leading only to exact re- 
lations among the CFs. According to the scheme of the 
composite operator method [ll|, it is possible to deter- 
mine these parameters by means of algebra constraints 



fixing the representation of the GF. By following this 
scheme, in Section IIIII we obtain extra equations clos- 
ing the set of relations and allowing for an exact and 
complete solution of the 11? BEG model. Section IIVI is 
devoted to the study of the finite temperature proper- 
ties. Finally, Sec. |V] is devoted to our conclusions and 
final remarks, while the appendix reports some relevant 
computational details. 



II. COMPOSITE FIELDS AND GREEN'S 
FUNCTIONS 

The Blume-Emery-Grifhths (BEG) model consists of 
a system with three states per spin. For first-nearest 
neighbor interaction the one-dimensional BEG model is 
described by the Hamiltonian 



iJ = - J^ 5(i)S'(« + 1) - A' ^ S-^ WS-^l^ + 1) 

i i 



(1) 



where the spin variable S{i) takes the values S{i) — 
— 1,0, 1. We use the Heisenberg picture: i = {i,t), where 
i stands for the lattice vector Ri. This model can be 
mapped into a fermionic model by means of the transfor- 
mation 



Si^) - [n{t) - 1], 



(2) 



where n{i) = X^o- ci(*)C(j(i) = c^{i)c{i) is the density 
number operator of a fermionic system; c{i) {c'{i)) is the 
annihilation (creation) operator of fermionic field in the 
spinor notation and satisfies canonical anti-commutation 
relations. Under the transformation ([2]), the Hamiltonian 
([1]) takes the form 

i i 

+ lwJ2ni^)D'^{^) + lwJ2Di^)n''{^) (3) 

i i 

i i 

where we have defined 

V = -{J + K), 

fl^-{2J + 2K-A-h), ^ ^ 

Hereafter, for a generic operator $(i) we shall use the 
following notation: $"(i) = [$(i + 1) + $(« - l)]/2. The 
Hamiltonian ([3]) is pertinent to an Hubbard model ex- 
tended to include intersite interactions, namely: charge- 
charge {V), charge-double occupancy (W) and double 
occupancy-double occupancy {—W) interactions. How- 
ever, Hamiltonian ^ does not exactly corresponds to 



the BEG Hamiltonian ([T]) since the mapping between S 
and n should take into account the four possible values 
of the particle density {n{i) — 0, n|(i) — 1 and ni{i) — 1, 
n{i) — 2). Letting the zero-state spin being degenerate, 
makes the Hamiltonians ([l} and ([3]) equivalent, provided 
one redefines the chemical potential fl and the on-site 
potential U as 



fL ^ 11 — p. ~ f3 "'^ln2, 
U ^ U = U -2/3-^\n2, 



(5) 



where /? = l/ksT. As a result, for a translationally 
invariant chain, one has 



2^ n(i)n"(i) + C/ 2^ D{i) + W^ n{i)D°'{i) 



-WY,D{i)D''{i) - Ai^ri(i) + Eq. 



(6) 



It is easy to see that the partition functions relative to 
the two models ^ and ([6]) are the same. Therefore, in 
the following we shall consider the fermionic model, de- 
scribed by the Hamiltonian ([S]), which is exactly equiv- 
alent to the BEG model. To solve this Hamiltonian we 
shall use the formalism of Green's functions and equa- 
tions of motion. As a first step, we show that there 
exists a complete set of eigenoperators and eigenvalues 
of H . To this purpose, one can introduce the Hubbard 
projection operators 



m = [1 - n{i)]c{i), 
rjii) — n{i)c{i). 

These fields satisfy the equations of motion 



(7) 



d 

i^v{i) = [-M + C^ + (2^^ + 2W)n"'{i) ~ WD°'{i)] rjii). 
at 

(8) 

By noting that the particle density n{i) and double oc- 
cupancy D{i) operators satisfy the following algebra 



n^(?) = n{i) + apD{i), 
DP{i)^D{i), 
nP{i)D{i) ^ 2D{i) + apD{i), 



for p>l, (9) 



where ap = 2^ — 2, it is easy to derive the following 
recursion rules 



-,a(i\]p 
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[D'^{t)Y 
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(10) 



The coefficients Am and Bm are rational numbers, satis- 



i(p) 



,(p) _ 



fying the sum rules J2m=i ^™ = ^ ^^'^ J2m=i ^^' = 1' 
whose explicit expressions are given in the appendix. 

On the basis of the equations of motion dS]) and by 
means of the recursion rules (fTO|) . it is easy to see that 
the composite multiplet operators 



V'^«H^) = 



V e,{i)[D-{i)]^ j 



i^^^^i) = 



J7(i)n°(i) 

')W[»°(»)]' 
V{i)D''{i) 



(11) 



are eigenoperators of H: 






(12) 



The energy matrices e^^' and £"' are matrices of rank 
7x7 and have the expressions 



e(«) = 
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where T = 2V + W. The energy levels are given by the 
eigenvalues of the energy matrices and are 



Ei^^ 



( -^ 

-/^+2y 
-Ai+vy+2y 

-/i+W/2+2y 

-^+vy+4y 
V-p+w/2+3y/ 



M''^ 



-M+t/+TV/2+y 

-Ai+(7+H'+2y 

-M+c/+2y 
-/i+c/+vy/2+2y 

-Ai+C/+W+4y 

V -p+c/+w+3y / 



(15) 

The knowledge of a complete set of eigenoperators and 
eigenvalues of the Hamiltonian allows for an exact ex- 
pression of the retarded Green's function 



G(^)(t - i') = d{t - i')(W^^H*, i), V^<^^^(*, i')}) 

i r+^ , (16) 



and, consequently, of the correlation function 
1 



+00 (17) 

(27r)y_oo ^ -^ 

In the above equations s = ^,r] and (• ■ ■ ) denotes the 
quantum-statistical average over the grand canonical en- 
semble. One finds 



G(^)(c.) = ^^ 



(s,ri) 



„^i Lu - £;^°^ + i6 



.('') 



p('*) 



where Tn = 1 + tanh(/?£'n /2) and the spectral density 
matrices a^^'"^' are computed by means of the formula 



^ fiv /^n/_^L nAJ At/ 



(19) 



A=l 



where $7^*^ is the 7x7 matrix whose columns are the 
eigenvectors of the matrix e^'^K The exphcit expressions 
of the spectral density matrices are given in the appendix. 
/('*■' is the normalization matrix defined as 

/(^) = ({^(^)(z),V>-^\*)}). (20) 

By means of the recurrence relations pO|) , all the matrix 
elements of I^^' can be expressed in terms of only the 
elements belonging to the first row. The calculations of 
the latter gives, for a homogeneous state, 



;^(fe-l) 



iR - -''-'' 



j(?) ^ §(k~b) _ Q(k~5) 



-in) 

l.k 

i-n) 

l.k 



^ik-1) 
Mk-5) 



where 

k(p) = ([n"'{i)]P), 

A(^) = lMz)K(z)f), 



(fc = l,...5), 
(A; = 6,7), 

(/c = l,...5), 

(fc = 6,7), 

S(P) = {[D'^{i)]P), 



(21) 



0<^P^ = -{n{^)[D'^{^)]P). 



(22) 

In conclusion, in this Section we have shown that the 
ID BEG model is exactly solvable. Exact expressions for 
the GF and CF have been obtained and are expressed in 
terms of a set of local correlation functions (P^ . which 
must be calculated in order to obtain quantitative re- 
sults. This problem will be considered in the next Sec- 
tion, where a self-consistent scheme, capable to compute 
the internal parameters, will be formulated. 



III. SELF-CONSISTENT EQUATIONS 

On the basis of the computational framework pro- 
vided in the previous Section, it is evident that the 
GF and the CF depend on the internal parameters: /i, 
k(p) and A(P) {p = 0, ...,4), S'^p^ and 0^p^ {p = 1,2). 
For a homogeneous state (i.e., translationally invariant: 
(n"(i)) = {n{i)) and {D{i)) = {D°'{i))), there are twelve 
parameters to be self-consistently computed in terms of 
the external parameters n, V, U, W and T. A first set 
of self-consistent equations is given by the algebra con- 
straints 

Ct(«)4(«) + Vl{i)v\{i) = 1 - "t(«)> 

ei(0el(0+»7i(0'7l(*) = l-"i(*), 

V1{^)v\{^) ="i(i) - D{i), 



(23) 



from which one gets the following self-consistent equa- 
tions 

^iS + ^S = ^'""'^ - ^'""'^ (fc = l,-5), (24) 
^S + C|^2 = '5^'"'^ - ^^'"'^ (fc = 6, 7) , 



where the CFs in the l.h.s. of Eq. (1241) can be computed 
by means of the formula 



, 7 



(25) 



ra=l 



Equations (|24|) provide one with eight self-consistent 
equations. To determine all the parameters one needs 
other four equations. These can be derived by means of 
the algebra constraints 



^t(i)„(i) = ^\i)D{i) = 0. 



(26) 



By exploiting these relations one can express the CFs 
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.(«) 
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Cll"' and Cll"' as 
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c[f(lx, + x, + lx! 



df) (ixi + -X2 + -XiX2 + -X2 

11 1^ i'2 ^2 4 ^ 



clf = cii^' (>i + Ix. + lx,x, + \xl + \xl 



'-'17 — '-'11 



1 



:X. 



1 



:Xi 



(27) 



The two parameters Xi and X2 are expressed in terms 
of the CFs C[l^\ C[l^^ and CJf^ as 



Xi = 



'-'12 

'-"11 



x.= 



C- 



(«) 



'-'11 



(28) 



Equations ([241), ^\ and (HSl) provide twelve self- 
consistent equations which will determine all the un- 
known internal parameters and therefore the various 
properties of the model. Once the parameters of the 
fermionic model are computed, by use of the mapping 
transformations (0]) and dS]) , it is straightforward to study 
the behavior of relevant properties of the BEG model. 
Details of the computations leading to Eq. (|77)) will be 
given elsewhere [I61] . 



IV. MAGNETIC AND THERMAL RESPONSES 

As an application of the general formulation provided 
in the previous sections, here we shall study the magnetic 
and thermal properties of the model, by restricting the 
analysis to the case ivT = and J < 0. In the following 
we set J = —1 and we consider only positive values of h, 
owing to the symmetry property of the model under the 
transformation h ^> —h. When A = the ground state 
is cither antiferromagnetic or ferromagnetic, depending 
on the value of the external field. As a consequence, the 
magnetization, defined as 



(5(z)) = (n(z))-l = l-2[c(f +C(?] 



(29) 



K=0 












:.::. 


h=2 




^^ 


-- 


— 


— ■ 








° 


T 




' 









(b) 














Q 


K=0 


















A=0 

A =0 25 


















^ 






h 


4 






















li=2 




/ K=0 





~ 


J=-1 






A^O.5 









7- 




' 



FIG. 1: (a) The magnetization m as a function of tiie external 
field hioT K = 0, J = ~1,T ^ 0.01 and for A = 0, 0.5 and 1. 
(b) The magnetization as a function of the temperature for 
iC = 0, J = -1, A = 0.5 and /i = 1, 2 and 3. 



FIG. 2: (a) The quadrupolar moment Q as a function of the 
external magnetic field h for K = 0, J = —1, T = 0.01 and for 
A = 0, 0.5 and 1. (b) The quadrupolar moment as a function 
of the temperature for K — 0, J = —1, A = 0.5 and h = 1, 2 
and 3. 



presents at T = two plateaus as a function of the exter- 
nal field h, with hg — 2 being the value of the saturated 
field. Turning on a positive single-ion anisotropy, three 
plateaus appear at rn = 0, m = 1/2 and m = \. This is 
in agreement with the criterion derived in Ref. [13] for 
the appearance of plateaus in spin chains in a uniform 
magnetic field. 

In Fig. [1^ we plot the magnetization as a function 
of the magnetic field at T = 0.01 for different values of 
A. Upon increasing the field, a nonzero magnetization 
begins at the critical value of the field he- this critical 
value decreases by increasing A. For h > he a m = 1/2 
plateau is observed until h reaches the saturated value 
hs, at which the third magnetization plateau at tti = 1 is 
observed. The width of the m = 1/2 plateau augments 
by increasing A in the range < A < 1 and becomes 
independent of A when A > 1. The critical field he and 
the saturated field hg satisfy, in the range < A < 1, 
the laws: 



2 -A, 



2 + A. 



Our findings are in goo d ag reement with Monte Carlo 
[14| and transfer matrix [I^ results. In Fig. [TJs we plot 
the magnetization as a function of the temperature for 
values of the magnetic field belonging to the three dif- 
ferent plateaus. For h — 1 {h < he) the magnetiza- 
tion is zero at T = 0; all spins are aligned (upward and 
downward) with the magnetic field, resulting in a pure 
AF state. When the temperature increases, the thermal 
fluctuations allow some of the downward spins to rotate 
and the magnetization increases up to T « 0.9 where it 
exhibits a maximum. Further increasing T, the thermal 
fluctuations enter in competition with the magnetic field 
and the magnetization decreases. For h = 2 (intermedi- 
ate phase he < h < hg) the magnetization is equal to 1/2 
at T = 0: half of the spins are parallel to h and half lie in 
the transverse plane. When the temperature increases, 
there is a slight increase of m (up to T ss 0.2), but soon 
the disorder induced by thermal fluctuation prevails and 
m decreases. For h = 3 {h > hs) at T = one finds 
TO = 1: all spins are parallel to the magnetic field and 



the system is in a pure ferromagnetic state. When the 
temperature increases, the long-range order is destroyed 
and the magnetization decreases. 

To further analyze the magnetic behavior, we have 
studied the quadrupolar moment Q, defined as 



Q={S^iz))^l 



2C[f. 



(30) 



At zero temperature also this quantity shows plateaus 
for A > 0. In Fig. [5^ the quadrupolar moment Q is 
plotted as a function of the external magnetic field, for 
K ^ 0, J = —1, T = 0.01 and for various values of A. 
Q takes the value 1/2 in the range he < h < hg, whereas 
is equal to 1 for all other values of h. The behavior of 
Q as a function of the temperature is shown in Fig. [2]d. 
For h < he and h > hg, the quadrupolar moment Q is 
maximum (Q — 1) at T = and decreases by increasing 
T. For he < h < hg, Q vanishes at zero temperature and 
increases augmenting T. 

The existence of the magnetic plateaus is endorsed 
by the peaks found in the magnetic susceptibility x — 
dm/dh. As evidenced in Fig. [3^, for A = one finds 
only one peak, whereas for A > there are two peaks 
appearing at he and hg, signalling a step-like behavior 
of the magnetization. When plotted as a function of the 
temperature, x shows a peak at low temperatures and 
then vanishes for T ^ for all values of the magnetic 
field but at hc^g, where, of course, it diverges due to the 
step encountered by the magnetization. In Fig. ^, we 
plot, as an example, the susceptibility in the neighbor- 
hood of hg. 

The specific heat is given by C = dE/dT where the in- 
ternal energy E can be computed as the thermal average 
of the Hamiltonian ([T|) for K — 



E = -J{S{i)Sii)°') + AQ~hm. 



(31) 



The specific heat exhibits a rich structure in correspon- 
dence of the critical values of the magnetic field. The 
behavior of the specific heat as a function of the tem- 
perature is shown in Figs. [3^-b in the neighborhood of 
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FIG. 3: (a) The susceptibility x as a, function of the external 
field for iC = 0, J = -1, T = 0.01 and for A = 0, 0.5 and 
1. (b) The susceptibility as a function of the temperature for 
K = Q, J = —1, A = 0.5 in the neighborhood of hs- 



FIG. 4: (a) The specific heat C as a function of the temper- 
ature for K = Q, J = —1, A = 0.5 and in the neighborhood 
of he- (b) The specific heat as a function of the temperature 
for K = G, J — —1, A — 0.5 and in the neighborhood of hs- 



hc^s- The possible excitations of the ground state are 
flipping of the spins parallel to h and/or of the ones per- 
pendicular to it. Far from the critical values, the specific 
heat presents only one peak at low temperatures. When 
h < he, the possible excitations are due only to the flip- 
ping of the longitudinal spins. In the neighborhood of he, 
a second peak appears since thermal fluctuations tends 
also to flip the spins perpendicular to the external fleld. 
Further increasing h, the specific heat presents only one 
peak until h k, hs, where again two peaks are present in 
a broader region with respect to h^. Away from h^, the 
specific heat shows only one peak. 



CONCLUDING REMARKS 



We have evidenced how the use of the Green's function 
and equations of motion formalism leads to the exact so- 
lution of the onc-dimcnsional BEG model. Our analysis 
allows for a comprehensive study of the model in the 
whole space of parameters K, J, A, h and T. Here, 
we have focused on the antiferromagnetic properties ex- 
hibited by the model and we have shown that, at zero 
temperature, the model exhibits three magnetic plateaus 
when A > 0. Furthermore, the specific heat shows a dou- 
ble peak structure in the neighborhood of the endpoints 
of the intermediate plateau. 
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APPENDIX A: SOME ANALYTICAL 
EXPRESSIONS 



by: 



The coefficients Am and Bm in Eq. (jTU]) are given 
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(Al) 



3 • 2P-1 



and 



B\ 



(P) _ o2-y 



B. 
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2 - 22-P. 



(A2) 



fj{s,m.) ^ sJ^'pC™) are the spectral density matrices de- 
fined in Eq. P^ : pC™) are matrices of rank 7x7: 

^1 k ^(1 0) 

rffe =(1 1/2 1/4 1/8 1/16 0) 
T-(3) 



l,k -(1111100) 

l,k 

l,k 

l,k 
-(7) 



= (1111111) 
=(11111 1/2 1/4) 



=(1 2 4 8 16 1 1) 



(A3) 



ik =(1 3/2 8/4 27/8 81/16 1/2 1/4). 
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and the S„i are given by: 
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m 



yis) _, 



2.7 — - 



-21 



-7I[%~7I[%+2I[%+3I[%-3I['>] 



'-3li% + lll[:>-12l[^l+4l[l 



Here we have reported only the first row of the spectral 
density matrices. All the other matrix elements can be 



expressed in terms of the first row by means of the recur- 
sion relation PH)) . 
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